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Abstract We study the effect of an imposed magnetic 
field on the motion of charged dust particles in mag- 
netically active regions of a protoplanetary disc. As- 
suming a power law structure for the vertical and the 
toroidal components of the magnetic field for the re- 
gions beyond magnetically dead region of the disc, the 
radial and the vertical velocities of the charged parti- 
cles, in the asymptotic case of small particles, are calcu- 
lated analytically. While grains with radii smaller than 
a critical radius significantly are affected by the mag- 
netic force, motion of the particles with larger radii is 
independent of the magnetic field. The critical radius 
depends on the magnetic geometry and the charge of 
the grains. Assuming that a grain particle has one ele- 
mentary charge and the physical properties of the disc 
correspond to a minimum-mass solar nebula, we show 
that only micron-sized grains are affected by the mag- 
netic force. Also, charge polarity determines direction 
of the radial velocity. For such small particles, both 
the radial and the vertical velocities increase due to the 
magnetic force. 
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1 Introduction 

Physical structures of an accretion disk are determined 
by various physical factors such as gravity of the central 
star, turbulence, magnetic fields. Understanding the 
time evolution of dust particles in circumstellar disc-like 
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structures around protostars and young stellar objects 
is crucial for explaining the formation of planetary sys- 
tems. Due to various complex physical processes that 
are playing significant roles in dust and gas discs, it is 
difficult to present a complete and self-consistent model 
for such systems. Because of the gas drag force, the mo- 
tion of the dust particles is different from that of gas 
(e.g., WeidenschiUing 1977). The gas drag force takes 
angular momentum from particles, resulting in their 
inward migration (Whipple 1972). Garaud, Barriere- 
Fouchct & Lin (2004) found that small particles move 
smoothly to the midplane, while large bodies oscillate 
about the midplane with decreasing amplitude. 

Takeuchi & Lin (2002) (here after TL) studied grains 
motion in optically thick disc, using the assumption 
that radiation pressure on the dust is negligible. They 
showed that as the gas accretion proceeds, small dust 
particles are left behind in the disc, and the dust to 
gas ratio increases. In a subsequent paper, Takeuchi 
& Lin (2003) investigated the outflow of dust particles 
on the surface layers of optically thick discs. Klahr 
& Lin (2001) studied the motion of grains in optically 
thin disc, where the dust experiences both gas drag 
and radiation pressure. In a similar analysis, Takeuchi 
& Artymowicz (2001) suggested that ring-like features 
which have been observed in some circumstellar de- 
bris discs (e.g. Schneider et al. 1999; Marsh et al. 
2005), arise naturally in discs passing through the tran- 
sitional phase from gas-dominated to dust-dominated 
discs. The dust particles migrate radially until they 
are in co-rotation with the gas. 

To our knowledge, most of the previous studies are 
assuming that the motion of particles in the gas and 
dust accretion discs exclusively caused by gravitation, 
gas drag force and radiation pressure of the central pro- 
tostar. However, coUisional dust grain electrification in 
the early solar nebula has long been discussed to be a 
cause for lightning, which could explain the melting of 
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mctcoritic chondrulc precursors. Poppc. Blum & Hcn- 
ning (2000) studied coUisional grain charging experi- 
mentally. They showed that coUisional grain charging is 
stronger than previously discussed with respect to pre- 
planetary grains and should be considered concerning 
the preplanetary dust aggregation and the formation 
of lightning in the solar nebula. It also has influences 
on other physical conditions in gas and dust discs, be- 
cause charged particles can interact with the magnetic 
field. Although the strength and direction of magnetic 
fields in accretion discs are rather unknown, the obser- 
vation of winds and outflows is explained by magnetic 
acceleration mechanisms that require the existence of a 
magnetic field with strength of ?» I0~* T at 1 AU solar 
distance (Konigl & Ruden 1993). 

The possible effects of magnetic fields on the struc- 
ture of accretion discs have been studied both nu- 
merically or analytically using simplified models (e.g. 
Lovelace, Romanova & Newman 1994; Bisnovatyi- 
Kogan & Lovelace 2000; Igumenshchev, Narayan & 
Abramowicz 2003). These studies, generally show that 
magnetic fields may have significant effects on the struc- 
ture of the discs. So, in gas and dust discs, one may ex- 
pect such type of behaviors due to the magnetic fields. 

However, the ionization structure of protoplanetary 
discs which determines the coupling between the mag- 
netic field and the grains is likely to be inhomogeneous 
and complex. The ionization fraction at different loca- 
tions, and even different times, are not easy to calcu- 
late (e.g. Ilgner & Nelson 2006). In fact, the idea of a 
magnetically dead-zone was first proposed by Gammie 
(1996). He argued that there is a poorly ionized re- 
gion where the growth of magnetorotational instability 
against ohmic dissipation can not be sustained. The 
dead zone typically stretching out to around 10 AU to 
15 AU (e.g. Reyes-Ruiz 2001; Matsumura & Pudritz 
2006). Wardle (2007) studied the ionization equilib- 
rium and magnetic diffusivity as a function of height 
from the disc midplane at radii 1 and 5 AU. He showed 
that the dead zone may not extend over the entire ver- 
tical dimension of the disc, because not only the ioniza- 
tion fraction but the gas density are strong functions of 
height above the midplane and so the magnetic coupling 
is likely to increase away from the midplane. Moreover, 
the ionization fraction itself is a strong function of the 
presence (and size) of the dust particles mixed with the 
gas and recent calculations have shown that once the 
grains have settled, the ionisation level is sufficient to 
couple the magnetic field to the gas at 1 AU, even at the 
midplane (Wardle 2007). But the disc becomes mag- 
netically active beyond the dead zone (e.g. Reyes-Ruiz 
2001). In this study, we are applying a simplified and 
phenomenological model for describing the gas compo- 



nent of a gas and dust disc for regions beyond the dead 
zone. 

We are assuming only the dust particles which are 
influenced by an imposed magnetic field. The proper- 
ties of any external magnetic fields threading a dust 
and gas disc are not well-known. However, it is very 
unlikely that the dipolar field of the central protostar 
is important at several AU from the protostar where 
the planets are believed to form. More probably such 
an external field comes from the the molecular cloud 
core out of which the gas and dust disc formed. 

These considerations show that magnetic forces on 
micron-sized particles in gas and dust disks will be 
worth further investigation. In fact, the structure of 
magnetic fields within protostellar discs may be studied 
via polarimetry provided that grains are aligned in re- 
spect to magnetic field within the discs. In this study, 
we assume that the charged particles are collisionless 
and indestructible spheres in a gaseous disc and, to 
keep the problem tractable, we also neglect the feed- 
back of the particles' drag on the motion of the gas. 
We mainly follow the approach of TL in analyzing the 
dynamics of grains, but in an imposed magnetic field. 
We show that a vertical imposed magnetic field affects 
the radial and the vertical motions of charged dust par- 
ticles. The basic equations of the model are described 
in section 2. In sections 3 we solve for individual parti- 
cles' velocities, which are charged and small with sizes 
smaller than the mean free path of the gas molecules. 
In the final section, we summarize our results. 

2 General Formulation 

We study the motion of a dust particle with mass m and 
charge g in a magnetized gas disc. One can simply show 
that the coUisional timescale is longer than the stopping 
time by at least an order of magnitude (Youdin & Shu 
2002). So, we can ignore the effect of collisions. For 
a charged dust particle, Newton's second law can be 
written as 

d^r 

-^grav Fg -|- Fii;^ag ^" -^electric i (1) 

where r is the radius vector of the particle. Also, Fgrav, 
Fg, F,nag and Fdcctric are the gravitational force of the 
central star, the gas drag, the magnetic force and the 
electric force, respectively. We neglect radiation pres- 
sure and Poynting-Robertson drag for simplicity. It is 
more convenient to write the equations of motion in 
cylindrical coordinates (r, 9) . 

As for the magnetic field geometry, we consider the 
toroidal and the vertical components of the magnetic 
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field. If a charged grain is placed in a rotating frame, 
it will invariably feel an electric field. The Keplerian 
velocity is denoted by wk = (GM*/r)^/^, where M, is 
the mass of the central star. We can imagine one of two 
approaches: (a) either state that the equation is written 
in a frame where electric field is zero, i.e. E + VK^g x 
B = in which case Lorentz Force term will become 
q{v — vkGo) X B; or (b) explicitly include an electric 
field term. We follow the first approach. In highly 
conducting plasma the electric field in the frame co- 
moving with the neutrals vanishes, so that the current 
density remains finite. But if we include the evolution 
of the magnetic field via the induction equation, this is 
not a good approximation in the presence of the charged 
dust grains. Because charged dust particles can easily 
decouple from the field lines. In order to avoid this 
difficulty, the magnetic field is treated as an imposed 
constraint in this study. 

Now we can write the radial and the azimuthal equa- 
tions of motion. 



dt 



r 



r 



{Vi.d - Wr,g) 



We now specify our ansatz for the magnetic field. 
The vertical component of the field is assumed to take 
a form 



r 



(6) 



where i?= 1 AU = 1.5x10^3 cmand^o = 10"-* T is the 
strength of the field at this radius. The normalization 
for the distance occurs relative to 1 AU. We will also 
consider more general values of n in order to explore the 
input parameter space. Such a power-law relation for 
the vertical component of the magnetic field has been 
studied theoretically in self-similar models (e.g., Shad- 
mehri & Khajenabi 2006; Wu & Lou 2006). For exam- 
ple, a distributed field with n = 5/4 has been sug- 
gested by Blandford & Payne (1982). The paraboloidal 
field model of Blandford (1976) corresponds to n = 1, 
and the models discussed by Ostriker (1997) corre- 
sponds to the range < n < 1. Also, the toroidal com- 
ponent of the magnetic field is (Lovelace, Romanova & 
Newman 1994) 



Be = t-Bb,s, 



(7) 



[{V9,d 

m 



dt 



{rve,d) 



-7jr{ve,d - ve,g) rv^.dBz, 



(2) 



(3) 



where Tg is the nondimensional stopping time. Also, 
Wr.d and vg^d are the r— and 6— components of the 
velocity of the dust particle, i.e. Vr^d = dr/dt = r and 
ve,d = rd9/dt = rO. Note that vg^g is the rotational 
velocity of the gas component. 

In order to write the z— component of the equation of 
motion, we assume the particles reach terminal velocity, 
at which the gravity, the gas drag and the Lorentz force 
balance with each other, 



-777fiK^ + qv^._dBg, 



(4) 



where fix ~ vyl/t. 

The nondimensional stopping time Tg can the be 
written as (TL), 



PgrvT '' 



(5) 



where wt represents the -^/S/tt times the mean thermal 
velocity c of the gas. Also, pd = 1.25 g cm'^ is the 
grain bulk density and pg is the gas density. In writing 
the above relation a spherical shape with radius s is 
assumed for the grain. 



where i?e,s = PeB^ and f3g is constant of order unity 
iPe < 0). ' 

The other physical variables of the gas component 
are 



T Z 

,(r, z) = Poi^T '^^P*^"2/J2 



cir) = co(-) 



(8) 



(9) 



(10) 



where the power-law indices p and q are generally neg- 
ative. Also, the radial and the rotational velocities of 
the gas are 



i),,g ^ -27ra(^)2(3p + 2q + 6 + 
K 



x(-i)^-t-V2AUyr-\ 



^]g = ^!K(l-77)l/^ 



5q + 9 z^, 
2 



(11) 



(12) 



where a is coefficient for the turbulent viscosity (ShakuraJ 
& Sunyaev 1973), and 77 is written as 



,/itr ,n , q + 3 z^ . 



(13) 
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For the overall disc properties, it is often sufficient to 
work with a disc model in which dust and gas are well 
mixed and coupled. In our analysis, we adopt the values 
correspond to a minimum-mass solar nebula: M* = 
IMq, po = 2.83 X fO-i" g cm-3^ ho = 3.33 x 10^^ 
AU, p = -2.25, q = -0.5 and a = 10"^. Having these 
values the gas mass inside 100 AU is 2.5 x I0~^ Mq. 



3 Analysis 

We can obtain approximate analytical solutions for the 
equations of the motion by applying some simplifying 
assumptions. From equations ([¥]) and ([7]), we obtain 



(14) 



where = {—q/'m)Bo is, in fact, the Larmor frequency 
which is defined for describing the motion of a collec- 
tion of moving charged particles in a region of constant 
magnetic field. 

We assume that rotational profiles of the gas and 
the dust components are close to Keplerian, i.e., vg_g ~ 
ve,d ~ wk- Thus, d{rve^d)/dt « Vr^dd{rvK)dr « 
Vr,dVK/2. Now, we can rewrite equation ^ as 



V0,d-ve,g = v,^dTs[—i — ) - -J 



(15) 



In equation ([2]), we can neglect the left hand side of 
the equation if (v^.d/vK)'^ 1. Thus, equation ^ can 
be written as 



+r—[{ve,d - wk)Sz - v^Bg] = 0. 



(16) 



Substituting from equation ()14p for d hito the above 
equation, an algebraic equation for Vi-^d is obtained and 
so. 



Vrd 



where 

I Ik r her r 



(17) 



(18) 




z/h. 
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Fig. 1 (top) Radial velocity Vr,d of dust particles of s = 
0.05 /^m, 0.1 /im, 1 fim and 10 /im at the radial distance 10 
AU from the central object. The distance from the midplane 
z is normalized by the disc scale hight kg. Dashed curves 
show profiles of the nonmagnetic solutions. The exponent of 
the magnetic field is n = 2. The electrical charge of a par- 
ticle is g = — e, where e is the elementary charge. For par- 
ticles with s = 0.1 /im the difference between the magnetic 
(solid curve) and the nonmagnetic (dashed curve) solutions 
is negligible. Profiles of the magnetic and the nonmagnetic 
solutions are the same for particles with larger radii (e.g. 
s = 10 /im). But for smaller particles (e.g. s = 0.1 and 0.05 
/im) the profiles significantly deviate from the dashed curve 
which shows nonmagnetic solutions for the same particles. 
{bottom) Radial velocity Vr^d of dust particles of s = 0.1 /tm 
for various exponent of the magnetic field n = 0.5, 1, 1.5 
and 2. Here, the dust particles are negatively charged with 
Q = -e. 



F, = 1 - 2-^(^)" + pli^fi^fi^f-. (19) 
Hk T Hk hg r 
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We can study profile of the radial velocity of dust 
particles using equation (jl7p . If we substitute this re- 
lation into equation (|14p . an analytical solution is ob- 
tained for the vertical velocity of the dust particles. 
Actually, for neutral dust particles, both parameters Fi 
and T2 become one and our analytical equations (|14p 
and ifTT]) reduce to the solutions of TL. The effect of 
the magnetic field on charged particles appears through 
Larmor frequency ft and the parameter f3g. However, 
in our proposed geometry of the magnetic field, we have 
an order of unity value for (3e. So, dependence of the 
solutions on this parameter is weak and in our sub- 
sequent analysis we assume Pe ~ ~1- But the Larmor 
frequency not only depends on the magnetic strength 
Bq, but mass and the electrical charge of a dust particle 
are important. We consider micron-sized grains where 
the nondimensional stopping time is much smaller than 
unity through most of the disc. Such small particles 
are well coupled to the gas component. Poppe, Blum 
& Henning (2000) demonstrated experimentally that a 
micron-sized grains may acquire 100 elementary charge. 
However, in cold protoplanetary discs, grains are either 
neutral or have 1 to 2 electronic charges. It is very un- 
likely planetary disc model where grains can carry 100 
charges near the mid-plane. It may be possible in the 
hot HII regions where temperature is very high, but 
in our analysis we consider 1 electrical charge. How- 
ever, this situation changes significantly as we consider 
higher vertical locations. For example, for r = 1 AU, 
the mean grain charge is — 14e above 3hg to 4/ig (Wardle 
2007). At larger radii, this transition may occur even 
closer to the midplane, given the drop in fluid density 
away from the central object. Since our solutions reach 
to z/hg = 3, it is plausible that at those locations the 
typical grain carries more than a single charge. 

We are interested in regions of the disc where mag- 
netically are active. As we discussed such regions must 
be beyond dead region which may extend up to 10 — 15 
AU. Although we consider the radial distance 10 AU in 
our Figures as starting radius of active region, we may 
simply consider larger radial distances in our analyti- 
cal solutions. Figure [1] shows radial velocity v^.d versus 
vertical distance z from the midplane of the disc for 
different dust sizes (top plot) and various exponent n 
of the magnetic field (bottom plot). In both plots the 
electrical charge of a dust particle is q = — e, where e 
is magnitude of the charge of the electron. In the top 
plot of Figure [U we fix the exponent n = 2. Also, solid 
curves correspond to the magnetic solutions and the 
dashed curves represent nonmagnetic solutions. For a 
particle with radius 1 /xm, the difference between mag- 
netic and the nonmagnetic solutions is negligible. Also, 
as the radius increases, the effect of the magnetic field 
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Fig. 2 The same as Figure (TJ but for positively charged 
particles with q = -|-e (solid curves) and the negatively 
charged particles with q = — e (dashed curves). 
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becomes weaker so that for a particle with radius 10 /im, 
magnetic and nonmagnetic solutions are indistinguish- 
able. But for particles with radii smaller than 1 ^m, 
the radial profile of the nonmagnetic solution is repre- 
sented by a dashed curve which is actually radial profile 
of the gas component as well. For such small particles, 
magnetic forces significantly modify profile of the radial 
velocity. In fact, the radial velocity of charged particles 
increases and as a particle becomes smaller, the radial 
velocity increases too. Note to the minus sign of the 
radial velocity. In other words, the radial migration of 
charged dust particles with radii smaller than critical 
Sc = 1 /im significantly increases due to the magnetic 
force. 

Bottom plot of Figure [T] shows radial velocity of a 
particle with radius 0.1 /im for various exponent n. The 
radial profile highly depends on the radial profile of the 
magnetic field. As the exponent n decreases, the radial 
migration of a charged grain increases. It implies that 
the critical radius Sc depends on the exponent n as well. 
While for n = 2 the critical radius Sc is around 1 /im, as 
the exponent n decreases, the critical radius increases 
which means particles with larger radii are affected by 
the magnetic force. 

Figure [2] shows the effect of the electrical charge po- 
larity on the radial velocity of the dust grains. Solid 
curves correspond to a positively charged grain with 
one elementary charge and the dashed curves are for a 
negatively charged grain. In this plot, the magnetic ex- 
ponent is n = 2. Although for grains with radii larger 
than Sc the effect of magnetic field is negligible, for 
grains smaller than this critical size the direction of the 
radial velocity depends on the polarity of grain's charge. 
Outward mass flux is carried primarily by the positively 
charged grains with radii smaller than the critical ra- 
dius Sc. Note that this result is based on the assumed 
magnetic field configuration. We considered a negative 
f30 according to Lovelace, Romanova & Newman (1994). 

Top plot of Figure El shows the vertical velocity v^^d 
of dust particles vs. z/hg for various sizes of grains 
and n = 2 and q = —e. Since the radial velocity sig- 
nificantly increases for grains with s < Sc due to the 
magnetic force, we expect larger vertical velocity com- 
paring to the nonmagnetic solutions according to equa- 
tion p4|) . Although both the radial and the vertical 
velocities increase due to the magnetic force for grains 
with radii smaller than Sc, sedimentation of the parti- 
cles can be determined based on the ratio of the time- 
scales of the sedimentation and the radial migration. 
The timescales of the sedimentation and the radial mi- 
gration are defined as tgcd ~ z/\vz,d\ and t^- w r/jwr.dl, 
respectively. Bottom plot of Figure [3] shows the ra- 
tio tr/tsed VS. distance from the midplane of the disc. 







\ \o.05 Jim 






0.1 \ 




1 




" n = 2 




\, . 1 \ 



3 


1 2 

z/ h 

a 








■' 0.1 ••r^'y''"'''^ 




/ 1 Jim ^^-^^--v. 






D 


1 2 

z/h 



Fig. 3 (top) Vertical velocity w^^d of dust particles versus 
the distance from the midplane of the disc for grains with 
radii s = 0.05 /im, 0.1 /im and 1 /im. The exponent of 
magnetic field is n = 2 and charge of a particle is q — — e. 
{bottom) Ratio of the time scales tr/tsed vs. distance from 
the midplane of the disc for the magnetic solutions (solid 
curves) and the nonmagnetic solutions (dashed curves). 



Dynamics of Dust Particles in Discs 



7 



Solid curves correspond to the magnetic solutions and 
the dashed curves are for the nonmagnetic solutions. 
Each curve is labeled by the radius of the grain. As the 
size of grain decreases, the ratio of timescales increases. 
When the effect of magnetic force becomes significant, 
the ratio of the timescales is larger than the nonmag- 
netic case. For example, grains with radius 0.05 /im 
may sediment before any large migration in the radial 
direction in the nonmagnetic case. But magnetic field 
force significantly enhances the ratio of the timescales 
to values larger than unity which implies large radial 
migration before sedimentation. Note that grains with 
radii 0.1, 1 and 10 fxra have positive radial velocities 
near to the midplane (see Figure [T]). However, as the 
distance from the midplane increases, the radial veloc- 
ity tends to zero at z = zq. Beyond this location the 
radial velocity is negative. On the other hand, the ra- 
dial migration timescale is proportional to inverse of 
the radial velocity. The dip at the curves occurs at 
the location where the radial velocity tends to zero (i.e. 
z ~ zo). But for grains with radius 0.05 /xm, the radial 
velocity is negative for all the distances from the mid- 
plane. So, there is not a dip in its corresponding curve 
of the ratio of the timescales. 



4 Discussion and Summary 

We presented a model for the dynamics of charged 
grains in a magnetized, protoplanetary disc. Properties 
of the gas component correspond to a minimum-mass 
solar nebula and the toroidal and the poloidal compo- 
nents of the magnetic field are assumed to be power-law. 
The vertical and the radial velocities of the charged dust 
particles depend on the geometry of the magnetic field, 
charge and mass of the grains. We assumed that each 
particle is charged one elementary charge. 

Our model illustrates the importance of charge of 
dust grains in their motions beyond the dead zone. Ob- 
servations of protoplanetary discs, at least in the opti- 
cal and infrared, mostly probe the upper disc layer. 
For charged grains with radii smaller than Sc, effect of 
the magnetic force on their motion is stronger in the 
upper layers comparing to the midplane region. Since 
magnetic force prevents sedimentation of such grains to 
the midplane of the disc, one may expect grains with 
radii larger than a critical radius in midplane region, 
unless smaller grains are neutral. In a disc resembling 
to the minimum-mass solar nebula, and assuming that 
all grains are charged, our model predicts that grains 
in midplane of the disc are mostly larger than Sc = 1 
fim. Dust sedimentation is known to affect the infrared 



spectra and the images of discs. In particular the far- 
infrared emission may be reduced by strong sedimen- 
tation. Although we did not investigate spectra of the 
disc, it is interesting to study possible modifications to 
the spectra due to the charge of the grains and the mag- 
netic force. However, in a more realistic model not only 
the coupling between the field lines and the charged 
grains should be considered but also the significant ef- 
fects of magnetic fields on the gas component. 

Takeuchi & Lin (2003) studied the outflow of dust 
particles in the upper layers of optically disc. They 
showed that grains in surface layers are moving out- 
ward due to the stellar radiation. The outward mass 
flux is carried primarily by particles of size 0.1 ^m ac- 
cording to their calculations. Although we neglected 
the stellar radiation in our model, we expect strong out- 
flows of sub-micron sized particles, if they are charged 
positively. In model of Takeuchi & Lin (2003), radi- 
ation pressure is the main physical factor in making 
surface outflows of grains. But in our model, it is the 
magnetic force which drives positively charged grains 
smaller than Sc to move outward. The actual optical 
depth along the radial direction in the midplane of the 
disk is huge, and thus, the radiation pressure from the 
central star onto the dust grains should be negligible. 
But neutral grains smaller than the critical radius in the 
midplane region can move outward once they charged 
positively. When they reach to the surface layers, not 
the magnetic force but the radiation pressure contribute 
to their outward motion. It will be interesting to study 
surface outflows under combined action of gas drag, the 
radiation pressure, and the magnetic force in order to 
understand their relative importance in making surface 
outflows. 

We can summarize our results: 

(i) Motion of particles with radii smaller than Sc are 
affected by the magnetic field, but larger particles move 
independent of their charge and the magnetic force. For 
input parameters corresponding to a minimum-mass so- 
lar nebula, the critical radius is around 1 /im and so, 
sub-micron sized particles are affected by the field lines. 

(ii) However, if the exponent of the magnetic field 
decreases or the charge of the grains increases, the crit- 
ical radius Sc increases. 

(iii) In our model, magnetic force increases the radial 
migration of charged particles comparing to when they 
are neutral. Also, vertical velocity increases due to the 
magnetic force. 

(iv) Although neutral dust grains tend to settle to- 
wards the midplane of the disc, magnetic force prevents 
charged grains with radii smaller than the critical ra- 
dius from settling down. 

(v) Direction of the radial motion is determined by 
the polarity of charged grains. 
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